Dynamics of bright matter- wave solitons in inhomogeneous cigar-type Bose-Einstein 

condensate 
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We discuss the possible observation of a new type of standing nonlinear atomic matter wave in 
the condensate: the nonlinear impurity mode. It is investigated dynamical effects of a bright soli- 
ton in Bose-Einstein-condensed (BEC) systems with local space variations of the two-body atomic 
scattering length. A rich dynamics is observed in the interaction between the soliton and an inho- 
mogeneity. Processes as trapping, reflection and transmission of the bright matter soliton due to 
the impurity are studied considering an analytical time-dependent variational approach and also by 
full numerical simulation. A condition is obtained for the collapse of the bright solitary wave in the 
quasi-one-dimensional BEC with attractive local inhomogeneity. 

PACS numbers: 03.75.-b, 42.81.Dp, ll.10.Lm, 42.50.Md 



I. INTRODUCTION 

Theoretical investigations of nonlinear collective exci- 
tations of matter waves, actually became a very interest- 
ing and relevant subject with the experimental observa- 
tions of Bose-Einstein condensation in vapors of alkali- 
metal atoms [1]. One of the interesting forms of localized 
waves of atomic matter are the solitons - moving sta- 
tionary nonlinear wave packets. Historically, it was first 
observed the dark solitons, represented by objects with 
nontrivial topological properties. They can exist in BEC 
systems for positive two-body scattering length (a s > 0) , 
corresponding to repulsive interactions between atoms. 
They can be presented as holes in the background of the 
condensates [2]. 

The observation of solitons in BEC systems with neg- 
ative two-body scattering length (a s < 0) - so called 
bright solitons, are more complicated from the experi- 
mental point of view. The difficulty is related to the 
instability of condensates in two-dimensions (2D) and 
three-dimensions (3D), when the number of atoms N ex- 
ceeds the critical limit N c , which is typically a number of 
the order of 1,500 atoms for Li 7 . This experimental limi- 
tation can be softened in case of quasi - one-dimensional 
(ID) geometry; i.e., for BEC in highly anisotropic cigar- 
type traps. We should observe that, for a true ID sys- 
tem, one does not expect the collapse of the system with 
increasing number of atoms [3,4]. However, it happens 
that a realistic ID limit is not a true ID system, with 
the density of particles still increasing due to the strong 
restoring forces in the perpendicular directions [5,6]. 

The observation of bright matter- wave solitons in BEC 
with attractive interactions has been recently reported in 
Ref. [7,8]. Theoretical models, explaining the observed 
phenomena, have been considered in Refs. [9,10]. It 



should be noted that in principle it is possible to observe 
the bright matter waves solitons in BEC with optical 
lattice. The possibility of changing the sign of the effec- 
tive dispersion in such lattices makes possible to generate 
bright solitons for repulsive condensates [11,12]. Thus, it 
represents an interesting possibility to control the dy- 
namics of bright matter-wave solitons. For the control, 
in this paper we suggest to use artificially induced inho- 
mogeneities, by considering a variation in the space dis- 
tribution of the atomic two-body scattering length. This 
variation can be achieved by using optical methods, as de- 
tuned laser field [13], or by means of Feshbach resonance 
applying an external magnetic field [14]. Mathematically, 
this leads to the appearance of a coordinate dependent 
coefficient in the nonlinear term of the Gross-Pitaevskii 
(GP) equation. We will analyze this problem for a cigar- 
type condensate using a full numerical solution of the 
GP formalism and also a time dependent variational ap- 
proach, which was successful in the description of BEC 
dynamics [15,5]. In spite of the well known difficulties of 
the variational approach it represents a good framework 
for a preliminary insight into the basic physical mecha- 
nism of the model. 

We consider the local space variation of the atomic 
scattering length as related to a nonlinear impurity term 
in the nonlinear Schrodinger equation (NLSE). Previ- 
ously, interaction of solitons with linear impurities have 
been performed in the sine-Gordon model [16,17] and in 
the ID NLSE [18-20]. It was shown that the interaction 
with the impurity mode leads to unexpected behaviors 
of the soliton. In particular, the soliton can be reflected 
by an attractive impurity This happens due to a reso- 
nant interaction between the soliton and an oscillating 
impurity mode. 
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For the description of the soliton dynamics, the varia- 
tional approach gives two coupled equations: for the soli- 
ton width and center-of-mass position. The oscillation 
frequencies of the width can be transferred resonantly to 
the oscillation frequencies of the center of mass; and, as 
a result, the trapped soliton can escape from the inho- 
mogeneity. We should also note that exists a solution 
representing a standing nonlinear atomic matter wave. 
The soliton, trapped by the nonlinear impurity, evolves 
to this solution. In this work we estimate the values of 
the parameters for the observation of such a new type of 
nonlinear standing atomic matter wave. 

The structure of the present work is as follows: In 
Sect. 2, we formulate the model for the quasi ID BEC 
with nonlinear inhomogeneity and derivate the equations 
for the soliton parameter using the time dependent varia- 
tional approach. In Sect. 3 we present the analysis of the 
fixed points, frequencies of oscillations for the width and 
center of mass. Section 4 contains the numerical model- 
ing of the system of ODE's and numerical simulations of 
the inhomogeneous GP equation, with a final discussion 
on the nonlinear impurity mode. 
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In this case, Z = Z(x 3 ,t) is normalized to the number of 
particles N. Next, we redefine Z and the variables, such 
that 
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From now, we consider a simplified notation for the 
functions, such that u = u(z,t), u t = du/dr, and 
u zz = d 2 u/dr 2 . Thus, we obtain the following ID NLSE: 



\u T = —u zz + er|it| 2 ?i — az 2 u, 



(7) 



where a is the signal of the two-body scattering length. 
The equation (7), for a = (uj^ — > 0) and a = — 1, has 
the solitonic solution 



V2 A sech [A(z - vt)} e 



vz ,n v 2 t 
~2 + 4~ 



(8) 



II. FORMULATION OF THE MODEL 

The mean field equation for a Bose-Einstein-condensed 
system, trapped by an harmonic potential, is given by the 
following GP equation: 



Airfi 2 a s , ,_ . l9 

+ -\W,*)\ 

m 



(1) 



where a s is the atomic scattering length, m is the mass 
of the atom, and the wave- function \1/ = ^>(f,t) is nor- 
malized to the number of particles N. In the present 
work, we assume a cylindrical highly anisotropic trapped 
potential, such that i.e. uj\ = ll>2 3> <^3- In this circum- 
stance, one can approximate the field as [5] 



*(r,t)=R(x u x 2 )Z(x 3 ,t). 



(2) 



where A is a constant and v is the soliton velocity. In 
the following, let us consider the interesting case of local 
space variation of a negative atomic scattering length, 
such that Eq. (7) is replaced by 



az 2 u- (1 + ef(z))\u\ 2 u. 



(9) 



e > (e < 0) refers to negative (positive) variation of the 
scattering length. The corresponding ID Hamiltonian 
energy is given by 



1 r c 

(H) = - 

n J-, 



dz 
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(10) 



where no is the normalization of u, related to the number 
of particles N and the scattering length a s : 



no = 47V|a s 



h 



(ii) 



R(xi,X2) = R satisfies the 2D harmonic oscillator equa- 
tion, that we assume is in the ground-state, normalized 
to one: 
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R = hmR, (3) 



R{xi,x 2 ) 



mwi 



(4) 



By substituting the Eqs.(2)-(4) in Eq.(l) and averaging 
over the transverse coordinates, 



We should observe that, in the realistic case one has a 
quasi-lD (cigar-like) trap. In Ref. [21], the authors have 
consider the formation and propagation of matter wave 
solitons, using a gas of 7 Li atoms, in a quasi-lD trap. The 
frequencies used in their trap are lu± = w\ = 2ttx 625 Hz, 
and Wi = W3 = 2irx 3.2 Hz, with the scattering length 
tuned to a s = — 3ao (ao is the Bohr radius). The ratio 
of such frequencies gives —a = (uj^/oji) 2 — 2.6 x 10~ 5 . 
In this shown in Ref. [6], the maximum crit- 

ical number no, ma x ~ 2.70 is a constant that does 
not depend on a s . The realistic maximum number of 
atoms N c , will be related to a s and the oscillator length: 
ANAas^mujx/Ti w 2.70. 
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To study the dynamics of the perturbed soliton we use 
the following trial function [22]: 



u = A sech 



(12) 



where A, a, (, <j), w, and b, are time-dependent variational 
parameters. In this case u is normalized to n = 2a A 2 . 
To derive the equations for the time-dependent param- 
eters of the soliton, we first obtain the averaged La- 
grangian 



L(t) 



J £(z,r)dz 



(13) 



with 



C(z,t) = -(u T u* -u*u) 



\u z \ 2 + az 2 \u\ 2 



(14) 



The equations for the soliton parameters are derived from 
the Lagrangian L, by using the corresponding Eulcr- 
Lagrangc equations. So, L is given by 



L = —no 
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We also obtain the coupled equations for a and (: 
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C TT = 4a( + e 



n dF 



(17) 



4a 2 d( ' 

When a is constant we have the well known description 
of the soliton center as the unit mass particle in an an- 
harmonic potential U(() = -It A 4 sech 4 (AC) [23]. From 
this point of view, for e < 0, a slowly moving soliton 
(v s < v c — 2-y/jej) can be trapped by the impurity. Con- 
sidering the impurity with e < 0, the soliton can pass 
through the impurity or reflect, depending on its veloc- 
ity. The numerical simulations of the system (17) and 
the GP equation (see section 4) shows that the dynamics 
is much more complicated. In particular we can observe 
the soliton reflecting from an attractive impurity even 
when v s < Ve- 
la order to have a more general formulation of the 
model, in the present section we have considered a non- 
zero external potential, parametrized by a. One could 
also explore the behavior of the soliton, by considering 
a more general time-dependent form of the external po- 
tential, as studied in Ref. [24]. However, in the present 
work our main motivation is the propagation of matter 
wave solitons, in a ID cigar-like trap [21], such that we 
will assume a — in the next sections. 



III. DYNAMICS OF BRIGHT SOLITONS UNDER 
TWO KIND OF INHOMOGENEITIES 



A. Point-like non-linear impurity 



Let us first consider an inhomogeneity given by a delta 
type (f(z) — 5(z)), and look for the fixed points of the 
system of Eqs. (17). In this case, 



F(a,C) =scch 4 (C/a), 

dF/dC = -(l/o) tanh (C/a) sech 4 ((/a) , 

dF jda = ((/a 2 ) tanh (C/a) sech 4 (C/a) . 



(18) 



The fixed point for the soliton center is given by C = 0. 
This corresponds to the case of an atomic matter soliton 
trapped by the local variation of the two-body scatter- 
ing length. In case the local variation corresponds to a 
positive scattering length (e < 0), we should observe the 
soliton being reflected by the inhomogeneity. Then, the 
stationary width a c can be defined by 



8 — 3e no 
2 n 



(19) 



Expanding the solution near a c we obtain the frequen- 
cies of small oscillations for the width a and for the 
center-of-mass C of the soliton, localized by the impurity. 
The square of such frequencies are, respectively, given by 



2 4 n 

uo = 



2 n 



7r 2 \8 — 3 no e, 



2 n 



^ \ 8 — 3 n e 



(20) 



In the variational approach for the soliton interacting 
with the impurity we have the interaction of the oscil- 
lating internal degree of freedom (the width) with the 
soliton center. As can be seen from Eqs. (20), the fre- 
quencies of the oscillations can match for a certain value 
of e, with energy transfer between the two modes. So, as 
a result of the reflection from attractive inhomogeneity 
in BEC, a depinning of the soliton can occur. 



B. Interface between two BEC media 



Now, we consider another interesting case of an inter- 
face between two media, such that at z — we have a 
sudden change in the two-body scattering length. The 
size of the inhomogeneity is given by e in Eq.(9), where 
f(z) = 0(z). In this case, we have 
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F(a, ()=a [2/3 + tanh ((/a) - ±tanh 3 ((/a)] , 
dF/da = F/a - (C/a)sech 4 ((/a) , (21) 
dF/dC = scch 4 (C/a) . 

And, from (17) with a = 0, we obtain the coupled equa- 
tions: 



16 

7T & 



(£ + 2)^^-6- 



tanlj | - 

a 



_ tanh J i + i sech 4 ± 
3 \a J \a 1 \a 

( TT = e ^seeh 4 (£) . 



(22) 
(23) 



There is no fixed point. At the interface, the value of the 
width is reduced, 



n (e + 2) : 



(24) 



and the frequency of oscillation of the pulse width is 

_ ng(e + 2) 2 

^ ~ 167T • 

For a constant value of a, from (23), we obtain 



(25) 



n 



2 rtp 
3a 



(26) 



When a = 4 /no the system (22) reduces to the single 
equation for ( describing the motion of the effective par- 
ticle, for optical beam crossing the interface of two non- 
linear Kerr media considered in [25,26]. If the velocity 
exceeds the critical value, the soliton pass trough the in- 
homogeneity. An interesting effect, predicted in Rcf. [26], 
can occur when the soliton cross the interface, namely 
the possibility of soliton splitting. The soliton is the so- 
lution in the first medium. In the second medium it can 
be considered as the initial wavepacket deviating from 
the solitonic solution, for this media. Applying the ap- 
proach developed in Ref. [27], such initial condition will 
decay on few solitons plus radiation [26] . The number of 
generated solitons is equal to 



n S oi 



i r 



\u \dz + i 



(27) 



where uo is the initial solution for the second medium, 
andX[...] stands for integer part of [...]. Thus, if we have 
a jump in the scattering length given by Aa s = a s2 — a 3 i, 
then the number of generated solitons in the second part 



\AWa s i + 1/2 



where 



of the BEC is equal to n so i = 

(1 + e) = a s2 /a s i is the ratio of the atomic scattering 
lengths. For example, for the ratio 9/4< a s2 /a s i <25/4 
(or 1.25< e <5.25) we obtain two solitons in the right- 
hand-side medium. To obtain n so i, we need e such that 
n so i{n so i - 1) < (e + 3/4) < n so i(n so i + 1). 



IV. NUMERICAL SIMULATIONS AND 
GENERAL DISCUSSION 



Our approach for pulses deviating from the exact soli- 
ton solution is interesting from the experimental point of 
view, considering the difficulty in producing exact soli- 
tonic solutions. Is of particular interest the non-trivial 
case of nonlinear Dirac-delta impurity (f(z) = S(z)), 
where we made detailed comparison between the vari- 
ational and full numerical solution of the GP equation. 
In Fig. 1, we are comparing the variation results with the 
numerical ones, for fixed point of the width, given by a 
(top frame); for the frequency of oscillations of the width, 
u> a (middle frame); and for the frequency of oscillations 
of the center-of-mass, uj^, trapped by the inhomogencity 
(bottom frame). By using the variational approach, we 
observe that the width goes to zero and the frequencies 
are singular when 



3 n 



(28) 



Here, it is interesting to observe that we have two criti- 
cal numbers that are related: one of the critical number 
no.max ~ 2.7 comes from the quasi-lD limit of a 3D cal- 
culation [6]; another, is the maximum amplitude of the 
delta-like impurity that we have just introduced, given 
by (28). Considering both restrictions, we have that the 
smaller value of e c is about one. 

The plots in Fig. 1 are valid for any value of no (where 
the maximum is about 8/3, according to [6]), because 
the width and the frequencies were rescaled, such that 
e' = (no/4)e (implying that e' c — 2/3), and 



, n ( 3 , 
a = —a = 1 — -e 
4 V 2 



no 



4/tt 



V(l-3e'/2) S 
no) (l-3e'/2) 2 ' 



As shown in Fig. 1, the variational results are supported 
by the full numerical calculation. The singularity occurs 
when the contributions coming from the inhomogeneity 
and nonlinearity are equal to the contribution from the 
quantum pressure, as seen from Eq. (17). When e > e c , 
occurs the collapse of solitary wave. So, we can observe 
the collapse of a ID soliton on the attractive nonlin- 
ear impurity. This possibility can be obtained following 
a dimensional analysis in the ID Hamiltonian given in 
Eq. (10). The behavior of the field at small widths is 
u <~ 1/L 1 / 2 . Taking into account that 8(z) ~ 1/L <~ |u| 2 , 
we can conclude that the contribution of the potential 
energy due to the impurity is ~ |u| 6 . For positive e, this 
term on the impurity exceeds the quantum pressure and 
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leads to the collapse of the soliton. In real situations, the 
collapse will be arrested on the final stage of the evolu- 
tion, when the soliton width become of the order of the 
inhomogcneity scale. Then, the delta-function approxi- 
mation for the impurity will break up. 




FIG. 1. The width a = o(no/4), the frequency of the width 
oscillations co' a = ui a (A/no) 2 , and the frequency of the soliton 
center oscillations u>£ = u^(4/no) 2 , versus the strength of the 
nonlinear delta- like impurity e = e(no/4). Solid line cor- 
responds to full numerical solution of the GP equation, and 
dotted line to the corresponding variational approach. All the 
quantities are in dimensionless units, as explained in the text. 




FIG. 2. Density profile evolution for a fixed value of the 
amplitude of the delta-like impurity, e = 0.3, in a projected 
3D plot. At z — 0, we observe the oscillation of the ampli- 
tude, starting from the normal one and going to the nonlinear 
localized one. Each line represents a fixed value of r. The plot 
is shown for a moderately small value of e, in order to avoid 
crossing of the lines. For larger values of e we can also ob- 
serve the emission of radiation, r and z are dimensionless 
quantities, as given in the text. 



In Fig. 2, we show numerical simulations of the wave 
profile. We note that, after strong emission of radiation, 
it evolves into the so-called nonlinear localized mode. The 
nonlinear localized mode represents an exact solution of 
GP equation with nonlinear impurity (9) and it is the 
nonlinear standing atomic matter wave. The solution is 
given by 



u (ni) = y/2asedi[a\z\ + f3] e 1 ' 



where 



= /3(e, a) = sign(e) In ^2|e|a + v 7 ! + 4e 2 a 2 ) 



1/2 



(29) 



(30) 



This solution can be obtained by using the solution of 
the homogeneous equation with the requirement of the 
field continuity at the inhomogcneity and satisfying the 
jump condition in the first derivative [28]. The normal- 
ization N( nl \ related to the umber of atoms for this so- 
lution is 



N (m) = 4a[l - e 7 ], 7 = 7 (e,a) 



7 = 



VI + 4e 2 a 2 - 1 
2e 2 a 



2 a 



VI + 4e 2 a 2 + 1 



(31) 



For small amplitude (or small impurity strength |e|) 
we obtain 



w 4a(l - e a). 



(32) 



At large amplitudes, we have -> 2/e for e > 0; and 

-> 8a for e < 0. Note that for e < we have a 
solution with two bumps structure for the nonlinear lo- 
calized mode. As shown in [28], this mode is unstable. 
Here, we have considered only the case e > 0. 

In order to verify the stability of the solutions, one 
can study the behavior of the second time derivative of 
the mean-square radius, as in Refs. [4,28,29]. To obtain 
the second time derivative of the mean-square radius, 
we use the Virial approach, with H = —d zz + V and 
V=-(l + eS{z))\u\ 2 : 

(z 2 ) TT = 4([H, zd z ]) = 8((-d zz )) - 4(zV z ), 



(Z 2 )rr = — I dz{S\u z \ 2 - 2\u\ 4 ) - —\u^. 

n J n 



(33) 



For the system to collapse we need (z 2 ) TT < 0; implying 
that 



e > 



2|u(0)| 



hwh 



|w| 4 )dz. 



(34) 



Using our solitonic ansatz, when a — > 0, we reach the 
critical limit, e c = 8/(3no) that was obtained before. 
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FIG. 3. Numerical simulations of the full GP equation, 
showing the dependence of no, related to the number of atoms 
N (top frame), and final velocity Vf (bottom frame), with re- 
spect to the initial velocity The results of both frames 
are shown for different values of e, as indicated inside the top 
frame. As shown, no = 4 was considered the initial value (for 
Vi = 0) of no. All the quantities are dimensionless. 



FIG. 4. For a fixed value of e = 0.4, it is shown the time 
evolution of the center-of-mass position, considering different 
values of the initial velocity, as given inside the frame. The 
results were obtained by using full numerical solution of the 
GP equation, considering no = 4. All the quantities are in 
dimensionless units. 



We have also investigated the dynamics of the matter 
soliton interacting with inhomogencity, studying differ- 
ent regimes of propagation for several values of e. In Fig. 
3, we present the results of numerical simulations for the 
final velocity (vf) versus the initial velocity (vi) of the 
soliton, considering different strengths e for the inhomo- 
geneity. In the present and next numerical approaches, 
we have considered uq = 4, for the general ID NLSE with 
nonlinear impurity, as one can easily rescale the obtained 
data to a value of no smaller than 2.7, in agreement with 
the quasi-lD results. 

As observed in Fig. 3, exists a region for the veloci- 
ties where the attractive nonlinear impurity reflects the 
soliton. In the model involving the constant width ap- 
proximation, this region corresponds to the trapped soli- 
ton. The numerical results show that always exists one 
window corresponding to the reflection of the soliton. By 
increasing e, this window is shifted to larger initial veloc- 
ities. From the variation of the number N, with respect 
to the initial velocity Vi (top frame of Fig. 3), one can 
also observe strong wave emissions by soliton, when e 
increases and tends to the critical value (sec also [28]). 
This picture reminds the picture of the collapse in 2D 
BEC with attractive interaction. We note that, by con- 
sidering the interaction of sine-Gordon kink with attrac- 
tive defect, many windows were found, corresponding to 
a resonance with local mode (see [16]). 

In order to compare with the results given in the lower 
frame of Fig. 3, we present in Fig. 4, for a fixed value 
of e = 0.4, full numerical calculation of the time evolu- 
tion of the ccntcr-of-mass position, considering different 
values of the initial velocity. It is clearly seen the repul- 



By considering numerical simulations of the variational 
equations (18) we found, qualitatively, the same behavior 
as the one observed in the full solution of the GP equa- 
tion (see Figs. 3 and 5). In Fig. 5, we show the behavior 
of the final velocity as a function of the initial velocity, 
for different values of e. In particular, the variational 
equations also show the existence of one window when 
the soliton is reflected by the impurity. When e increases 
the window moves in the right direction, reducing the 
width. 




-1 1 1 ' 1 1 1 1 

12 3 

v i 

FIG. 5. Final velocity Vf versus the initial velocity Vi, for 
different values of e, as indicated inside the frame, using 
no = 4. The results were obtained from numerical simula- 
tions of the variational equations given in Eq. (17). All the 
quantities are in dimensionless units. 
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In distinction, on the GP equation, the variational 
equations shows a more complicated (probably chaotic) 
dynamics near the points where the regime of reflection 
starts or finishes. In this region, we also can see the 
more rare events with the transmission of soliton through 
the impurity. This observation resembles the phenomena 
observed at the interaction of the sine-Gordon equation 
kink with a local defect. The system of ODE's (two mode 
model) has a similar structure as our equation (18), show- 
ing chaotic behavior, leading always to a finite time for 
the period of the soliton trapping. This phenomenon is 
due to stochastic instabilities inherent for this dynamical 
system. The reason for this phenomenon is that the fi- 
nite dimensional system, like the one given by Eq. (18), 
cannot take into account the soliton radiation, that inter- 
acts with the defect. The effect of the radiation leads to 
the appearance of the damping in Eq. (18), that changes 
the long-time behavior of the system. In particular, the 
radiative damping can lead to the long-time regular dy- 
namics [17]. 

For a more detailed investigation of the interaction be- 
tween matter wave soliton with the nonlinear impurity 
mode, excited on the inhomogeneity, we need to develop 
collective coordinate approach like the one considered in 
Ref. [18]. This will be considered in future. 



V. CONCLUSION 



atoms) exceeds a certain critical value. This effect in 
quasi-lD BEC resembles the phenomena that occurs in 
2D BEC. Using the time-dependent variational approach 
we have described successfully both phenomena. 

For the case of a sudden variation of the two-body 
scattering length, represented by a nonlinear jump inho- 
mogeneity, using analogy with a nonlinear optical prob- 
lem [26], we presented the condition for the multiple 
bright matter soliton generation. 

Finally, we would like to emphasize that the problem 
we have studied in the present work has important ap- 
plication for the control of parameters of bright atomic 
matter solitons and for the generation of solitons in quasi 
ID BEC. 
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In this work, we have investigated the dynamics of 
bright matter wave soliton, in BEC systems with cigar 
type geometry and attractive interactions. The inhomo- 
geneities can appeared in BEC due to the existence of 
regions in space with different values of the two-body 
atomic scattering length o s . These variations can be 
achieved using, for example, the Feschbach resonances. 
Two kind of inhomogeneities in the spatial distribution 
of a s have been studied: local point-like and jump type. 
The first type has been modeled by a Dirac-delta func- 
tion, that will result in a modulation of the nonlinear 
term in the GP equation, corresponding to the so called 
nonlinear impurity in the nonlinear Schrodinger equa- 
tion. The second type corresponds to a sudden variation 
of the two-body scattering length, that affects the ampli- 
tude of the nonlinear term in GP equation. It correspond 
to the case that the BEC system is divided in two parts 
with different values of a s (see discussion of recent exper- 
iment in Ref. [13]). 

The present investigation of the local variation in 
space of the atomic scattering length shows that differ- 
ent regimes of the soliton interaction with the nonlinear 
impurity are possible. It is observed trapping, reflection 
and transmission regimes. The most interesting effect is 
the reflection of the atomic soliton by the attractive non- 
linear impurity. We have also verified the occurrence of 
collapse of the soliton on the attractive impurity, when 
the strength of the impurity (or the initial number of 
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